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We employed the method of virial expansion in order to compute the retarded density correlation 
function (generalized diffusion propagator) in the critical random matrix ensemble in the limit of 
strong multifractality. We found that the long-range nature of the Hamiltonian is a common root 
of both multifractality and Levy flights which show up in the power-law intermediate- and long- 
distance behavior, respectively, of the density correlation function. We review certain models of 
classical random walks on fractals and show the similarity of the density correlation function in 
them to that for the quantum problem described by the random critical long-range Hamiltonians. 



I. INTRODUCTION 

Multifractality of critical wave functions^ has recently 
emerged as a subject of not only fundamental theoretical 
interestP but also of experimental activities at a fron- 
tier of technological achievement^^HH jn surface imaging, 
sound propagation and driven cold atom systems. It 
could be also an important issue in describing the su- 
perconductor (superfluid) to insulator transitions in con- 
densed matter and systems of cold atomJ^, Kondo ef- 
fect in strongly disordered metald^ and other interaction 
and non-linear phenomena in systems with strong dis- 
order. Multifractal correlations is the property of wave 
functions near the Anderson localization transition in 3D 
systems with short-range Hamiltonians as well as in sys- 
tems with long-range |-ffr,r'| ^ \r — r'\~'^ hopping^ d 
is the space dimension. In the latter systems critical cor- 
relations of wave functions may be realized in a broad 
range of parameters. 

Usually multifractality is studied by considering the mo- 
ments of the local density of states at a given energy e 
or the correlation functions of such moments at different 
energies. However, it may show up also in the density 
correlation function: 



(£-£„ + iO) (e - - w - iO) 



(1) 



In Eq. ([T]) we denote by 4'n{r) the exact n-th wave func- 
tion at a site r which corresponds to the energy e„ and e 
is the fixed energy (which we set zero throughout this pa- 
per). The symbol (...) stands for the ensemble averaging. 
The Fourier transform 2?r,r' (t) of this correlation function 
("generalized diffusion propagator") describes spreading 
of the wave packet in space and time and is a direct gen- 
eralization of the pure classical object: the probability 
density to find a random walker at a point r at a time 
t after the beginning of random walk from the origin, 
r' — 0. The possibility of both classical and quantum for- 
mulation of the density correlation function makes this 
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FIG. 1. (color online) Generalized diffusion propagator 
(density correlation function) for critical eigenstates of long- 
range Hamiltonians Eq. ( 12 I with strong multifractality 
d2 <^ d. The multifractality is encoded in the intermediate- 
distance (1 < r/^o < {Eot)i) dependence .^^^hile the 
Levy flights exhibit themselves in the long-distance {r/£o 3> 
{Eot)i) tails r"^''; Iq, Eq are the lower length scale and 
the upper energy scale of multifractalitjEH, respectively. The 
overall functional form of the density correlation function 
obeys the scaling law r"'' S{r'^ /t). 



object convenient for quantum-to-classical comparison. 

The descender of the generalized diffusion propagator is 
the survival probability or the return probability P{t) 



(2) 



where po is the averaged DoS. 

This is the probability for a wave packet to survive for 
a time t at a point where it was created at t — or the 
probability for a random walker to find himself in the 
origin at a time t. For extended states P{t) decreases 
with time, albeit differently for the truly extended states 



2 



and the critical multifractal stateJ^^IHl 



requires 



P{t) cx 



r 



extended states 
critical states 



(3) 



In Eq.([3]), d2 < d is the non-trivial critical exponent, 
the correlation fractal dimension, that determines inverse 
participation ratio in a critical system of the size L: 



l2 = {j2\^ir)n^L 



-d2 



(4) 



A general form of the density correlation function, 
V^KjUj), was suggested by C halker in the momentum- 
tioiPEH. 



frequency representatior 



In the region of critical 



states it can be written as follows: 



27r/9o 



i-iuj + o) [fi\K\L^,^) + iy 



where po ^ N ^Y.n=o(^i'^ ^ en))|e=o is the averaged 
density of states at e = 0. 

The key element of this prediction is that at criticality 
(in the classical Dyson symmetry classes) there is only 
one scale cx a;3 which separates small- and large- 

momentum behavior, both of them being different limits 
of one and the same analytic function /(cc). The com- 
bination L-ii^ (X {—iuj + 0)~'^ in Eq.([5]) is dictated by 
the retarded character of the correlation function which 
should be a regular function in the upper half-plane of 
complex uj. 

Note that the original idea of Chalker was that of recon- 
ciling the one-parameter scaling of Gang-of-FouJi^ with 
the multifractal behavior. The one-parameter scaling is 
encoded in Eq.([5]): the counterpart of 'D{k,uj) in the 
coordinate-time domain obeys the scaling law 



(5) 



r''A.,o{t) = SirVt), 



(6) 



where S{x) is the scaling function. 

According to Chalker 's conjecture, the mult ifr act ality ex- 
hibits itself in the asymptotic form of the function f{x) 
at \x\ > 1: 



> 1. 



(7) 



Here c?2 is the same correlation fractal dimension defined 
by Eq. Q which appears also in Eq. ([3| as well as in 
the local density of states correlation function Ref.'^. It 
determines the multifractal regime of the density correla- 
tion function (see Fig. 1) where 



{t > r'^) . 



(8) 



At small I a; I ^ 1, the function f{x) is small as the prob- 
ability conservation 



J2'Dr,r'it) ^2T:poeit) 



(9) 



fix = 0)=0. 



(10) 



Assuming that a; = is the regular point of the analytic 
function f{x), one can expand 



f{x) — ClX + C2X^ + CsX^ + 



{x = \k\ L- 



(11) 



Note, however, that the regular expansion in x, Eq.(|ll[), 
does not mean analyticity in the momentum domain. For 
instance the linear in x term contains non-analytic term 
\k\. Non-analytic in k terms arise from any odd power of 
X. These terms imply power-law tails of the density corre- 
lation function in the coo rdin ate domain. Namely, if the 
leading odd term in Eq. (11) is x'^™~^, the density cor- 
relation function acquires a power-law tail oc 7-~(2™-i+<i) 
at large r. By analogy with classical random walks, we 
will refer to such tails as the Levy flights. 
The simplest class of critical systems, where such non- 
analytic terms are present, includes the systems with 
long-range resonance hopping where the hopping terms 
of the Hamiltonian decrease with distance as: 



(12) 



However, the situation depends a lot on the dimension- 
ality d. As we will show below, the Levy flights in such 
systems are characterized by a power-law tail: 



Vr.a{t) cx r 



-d 



t 



(t < r^) 



(13) 



For 0? = 3, the second moment (r^) does exist so that 
we have in Eq. (11) that ci = and C2 is finite. The 
non- analytic in k terms appear starting from x^, which, 
however, is sub-leading at small x. The role of non- 
analyticity is greatly enhanced at d = 1 where ci ^ 0. 
Here the leading term in Eq. (11) is non-analytic in k. 
The case d = 2 is speciaP^'. In this case even criticality 
is not guaranteed by Eq. ( 12 ) and the answer may de- 



pend on the particular symmetry class and details of the 
Hamiltonian. 

In this paper, we concentrate on the simplest system 
which belongs to the class of long-range random Hamil- 
tonians, Eq. (12), with d — 1. This is the ensemble of 



Hermitian critical power-law banded random matrices^ 
(CPLBRM) where the multifractality can be continu- 
ously tuned from weak to strong by changing only one 
parameter < b < oo: 



{Hnm) — 0, {\Hnm\^) 



n^m, (14) 



where /3 = 1, 2, 4 is the Dyson symmetry number. 
We consider the limit of strong multifractality h ^ 1 
and demonstrate the validity of Eqs. ([s]), ([t]), (|8|, ( [iT| ) 
and ( 13 ) which are the key relationships for the descrip- 



tion of Levy flights and multifractality. Furthermore, we 
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will find the function f{x) explicitly for the Gaussian Or- 
thogonal (GOE) and Unitary (GUE) ensembles thereby 
unifying both concepts of multifractality and Levy flights 
in the one single analytic function. 

One of the purposes of this paper is to compare the 
density correlations function for random critical states 
of a long-range Hamiltonian with that for the classi- 
cal random walks on a fractal. For a classical random 
walk with a constant step length on a lattice the den- 
sity correlation function decreases exponentially at large 

distances Vryit) oc e~2Dt and thus all moments (r ) 
are well defined. A peculiar class of random walks called 
"Levy flights" is realized in the models with variable step 
length. If the probability to have a long step is heav- 
ily tailed, the probability distribution of random walks 
may have power-law tails at large distances and thus the 
moments (r^™) may diverge^^. The similar power-law 
tails may appe ar in classical random walk on hierarchical 
manifold4iSll20|^ where fractality is build in their geome- 
try. 

We will demonstrate that certain models of classical ran- 
dom walk on a fractal result in the same regimes of the 
density correlation function as for the quantum critical 
long-range Hamiltonians of Eq. (12). Namely, the two 
principle regimes Eqs. ([8|, (13) described in Fig. 1 have 
their counterparts for the classical diffusion, with the 
correlation fractal dimension d2 being replaced by the 
Hausdorff dimension dh and the dimensionality of space 
d being replaced by the walk dimension. Moreover, we 
will see that the critical exponent 2/d oi sub-diffusion 
(r^) oc i^/"* for critical random Hamiltonians in d > 2 di- 
mensions appears to have the same relationship with the 
exponent —d of Levy flights r'^ 2?o,r(0 = S{r^ t) oc r~'^ for 
the quantum problem with long-range random Hamilto- 
nians, Eq. (|l2l), and for the certain classical random walks 
on fractals 

The paper is organized as follows. In Section II we review 
the method of virial expansion in the number of resonant 
levels and give general expressions for the density correla- 
tion function (effective diffuson) for any almost diagonal 
Gaussian Random Matrix Theory (RMT). In Section HI, 
we apply these expressions to compute the density corre- 
lation function for the specific CPLBRM in the coordi- 
nate/time representation and uncover multifractality and 
Levy flights as consequences of the long-range nature of 
the Hamiltonian, Eq. ( 14 ). In Section IV, we compute the 



density correlation function in the coordinate/frequency 
representation. Sections V and VI are devoted to ver- 
ification of the key relationships Eqs. ([5|, ([7| and (11) 
due to multifractality and Levy flights as seen in the 
momentum-frequency representation of the density cor- 
relation function. We explain a perturbative derivation 
of Levy flights for strong multifractality and the Wigner- 
Dyson limit of the CPLBRM in Sections VII and VIII, 
respectively. In Section IX, we present numerical results 
for the critical RMT with weak and strong multifractal- 
ity. In Section X, we discuss the classical random walks 
on fractals and show that in certain models they have the 



same phenomenology as their quantum counterpart. In 
Section XI we summarize the main results of the paper 
and discuss their applicability to the Anderson transition 
in 3D systems. 



II. VIRIAL EXPANSION METHOD 

In order to check the validity of Eqs. ([5]), ([t]), 
and (pl for the ensemble of CPLBRM, Eq. 




strong multifractality, 6 ^ 1, we exploit the meth od of 
virial expansion in the number of resonant leveld^HHl^ 
This method is a certain re-summation of the locator 
expansioi^^ in the hopping matrix element hnm'- 



1 



2 (n — ?T!,)2 



(15) 



The summation is organized so that any correlation func- 
tion C{K,io) is represented as a series in 6^1: 



where the length scale is itself 6-dependent: 



2\/26 

w/3 ' 



(16) 



(17) 



and 77 is specific to a particular correlation function. 
The advantage of this representation is that Eq. ( 16 ) is a 
functional series, each function C^™^\kL^) adding details 
of correlations which emerge from resonant interaction 
of m states. These states have energy levels £m being 
anomalously close to each other within the interval of or- 
der 5^1. At small 6, such "multiple collision" of levels 
has small probability. Therefore, even the few first terms 
in the series Eq. ( 16 ) give a very good approximation of 



the correlation function. 
Note that the density of states pe ~ \i 2-n 
responding to CPLBRM Eq. ( [l4| ) has negligible vari- 
ation at a scale E ^ b <^ 1^24J. We will neglect 
such variations throughout the paper and approximate 
p{E) ss p(0) = po- With such an accuracy we obtain for 
the retarded density correlation function: 



27r/9o 



(18) 



In Eq. ( 18 ) we denote 



2npoV^^;^\oj)= I e''-'dtV'^;V{t)dt. (19) 
Jo 

The density correlation function in the time domain is 
given by: 

-Dry it) ^ (2^/90 5ry + Vfl, (t) + Vfl, (t)) 0{t) . (20) 



it^t o 7^(2,3) / 
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To simplify formulae, we assume below t > and skip 

o{t). 

Calculating leading terms of the virial expansion, we ob- 
tain the following expressions for /? = 1 (GOE): 

it) = 27rhl,, t e-^^^y'^' /o {{Kytf) ■ (21) 

and for /3 = 2 (GUE): 



-ih.tf 



hr^r' Erf(/lr,r'^) 



(22) 

where hry is defined in Eq. (15); . .) is the modified 
Bessel function of the first kind. Answers for subleading 
terms of the virial expansion are more lengthy and we 
will publish them elsewhereF^. 



Equations ( 18 )-( 22 ) are valid for any Gaussian RMT with 
an arbitrary variance of the hopping matrix elements 
that is parametrically smaller than the variance of the 
diagonal ones [almost diagonal Gaussian RMT). 

Remarkably, Vi^^f} (t) at r = r' is expressed in terms of 
'Drjtkt) as follows: 

2^^"^ (0 = - E ^t' (^) , = 2, 3. (23) 

r^r' 

This is a consequence of the particle conservation which 
requires 

Y,A,r'it)=27Tpo. (24) 
r 

As this relation is already fulfilled by the first term in 



Eq. (181, all other terms proportional to V^.^; (to > 2) 
shouli 



d obey Eq. ( 23 ) 



The small-distance behavior Vryit) oc \r 



is also 



remarkable. Applying Eq. (231 one obtains for the sur- 
vival probability Eq. ([2]): 



P(t) = 1- ^ ln(bt) w (bty 
d 



(27) 



(see Ref.JT^ for more details) where for the CPLBRM 
one findpl 




GOE 
GUE 



(6 « 1) 



(28) 



Thus Eq. (27) shows that the asymptotic behavior \r 



is a signature of (strong) multifractality with the 



correlation dimension given by Eq. ( 28 ) 



We conclude this section by emphasizing once again that 
the multifractality and Levy flights for the long-range 
Hamiltonians Eq. ( 12 ) have the same root and cannot 



exist one without the other. 



IV. DENSITY CORRELATION FUNCTION IN 
THE FREQUENCY DOMAIN 



One can invert Eq. ( 19 ) and obtain 



where 



z^{\r-r'\/L^Y o^u\r-r'\ 



(29) 



(30) 



c? = 1 for CPLBRM, is given by Eq. ( 17) and expres- 
sions for read: 

« sgn(z) (/o(z2) - hiz"))] , (GOE, /3 = 1) ; (31) 



III. DENSITY CORRELATION FUNCTION IN 
THE TIME DOMAIN AND RETURN 
PROBABILITY 



An inspection of Eqs. (^2lh, (22) shows that for the 



CPLBRM with & ^ 1 the asymptotic behavior of the 
effective diffuson V^y it) = ^, (t) is given by: 



r V2b 



\t I 2 



GOE 

77, GUE 



\r ~ r'\ <^ht . 



(25) 

|r-r'|>6t. (26) 



In a mor e ge neral case of long-range Hamiltonians obey- 
ing Eq. (12) one should replace |r — r 



r - r'\'^ to 



obtain the large-distance tail Eq. ( 13 ) which is the hall- 
mark of the Levy flights. 



/ 

l + i \ Erfi(z) 



(32) 



/^(z^ + i) 



, (GUE,/3 = 2) 



Kj{...) is the modified Bessel function of the second 



kind. Functions ^'^^\z) in Eqs. (31 1, (32 1 are qualita- 
tively very similar for the orthogonal and the unitary 
ensemble (see Figs. |2|3l ). Their real parts are even func- 
tions of z which exponentially decrease like \z\^ e~'^'^ 
at l^l 3> 1. The imaginary parts are odd functions of z 
which behave like z~^ at large \z\. Correspondingly, the 
density correlation function in the coordinate-frequency 
representation behaves at \z\ ^ 1 (large distances or fre- 
quencies) like: 



Im[I?r^r' (w)] ^ sign(z) \z 
Re[Pr^r'(w)] ^ 



(33) 
(34) 
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Real part 




FIG. 2. (color online) Real part of the functions f7<°°'^'(2) 
in Eq. ml (blue) and n'-^^^\z) in Eq. ^ (red). They 



(2) 

describe the cj-dependence of the corrections ^, to the re- 
tarded density correlation function at a fixed r ^ r' . Note 
that these corrections are caused by the extended nature of 
wave functions. 



Imaginary part 




FIG. 3. (color onl ine) Imaginary part of the functions 
q(GOE)(-^-) in Eq. ^ (blue) and ^'^'^^'^\z) in Eq. ^ (red). 



Note that the correlation function 'D,.^r'{'^) for jr — r'| ^ 
Lijj has a power-law tail. This is another manifestation 
of Levy flights in the r — uj domain. 



FREQUENCY-MOMENTUM 
REPRESENTATION. 



Now we compute the retarded density correlation func- 
tion 'D{k,^u)) in the momentum-frequency representa- 
tion. This allows us to make contact with Eq. ^ to 
check Chalker's conjecture on the effect of multifractal- 
ity Eq. ([7| and to detect the non-analytic in n terms in 
Eq. ( [TT| ), which is a signature of Levy flights. 
To this end we define: 



(35) 



where N x N is the size of the matrix. We assume this 
size to be larger than all the relevant length scales of the 
system. If in addition k <C 1, one can switch from the 
sum to integral in Eq. ( 35 ) to implement the continuous 
approximation which is valid at \Lt^ \ ^ 1. 
Before proceeding with calculations, let us note that 
Eq. ([23l) implies: 



V^^^K = 0,u})= 0. 



(36) 



Thus one can do the Fourier transform of 'D^^,{uj) at 
|r — r'l 7^ and then subtract the k = value of the 
result obtained. Thus, we obtain from Eqs. (31), (32): 



-iuj + 



+ \L^\G{fi), 



(37) 



where n = k . One can easily see that Eq. ( 37 1 is 



consistent with Eq. ([5| in the limit 61og(fc) ^ 1, with 
f{x) — j(-\/2/7r/3po) ^G'(fc). The retarded character of 
the correlation function is encoded in real and imaginary 
parts of G{k) which obey the Kramers-Kronig relation- 
ship. The expressions for them read 



GOE, /3 = 1: 



Re[G(«o^=)(S)]-^|«:|e 



16 In 



32 



Im[G(°°'^)(K)] = -sign(cj)— x 
4pFq({l,l,3},{^,2,[ 



(38) 



pFq({l,l},^i,^ 



2 \ 1 



(39) 



GUE, I3^2: 
Re[G(«ui=)(;i) 



TT 

Y 



K e * 



(40) 



J _ 

Im[G(°UE)(-)] ^ _sign(w)^ [2V^e~T ^Erfi 



s2pFq({l,l},<i^,2 



2 \ 1 



(41) 



Here 7 is the Euler constant and pFq(...) is the hyper- 
geometric function. We discuss these results in the next 
section. 



VI. COMPARISON WITH THE CHALKER'S 
CONJECTURE. 

The asymptotic behavior of the functions G'^^^(k;) and 
is very similar. At small n <^1\ 



r,r' 



Mcm-l^Sl^; III 



(42) 



6 



Im[G(K)] = 



(43) 



8 



This is consistent with Eq. ( |ll[ ) (at d = 1) where all Ci ~ 
b. In particular, the non-analytic behavior cx \k\ caused 
by the Levy- flights is confirmed by Eq. (|42|. Eq. (43) 
corresponds to corrections oc C2 x"^ in Eq 
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FIG. 4. (color online) Real part of G(fc) for GOE (blue) and 
GUE (r ed). The dashed lines is the asymptotic behavior given 
by Eq. ([Tl]) (red) and Eq. ^ (black) at = O.Of and c = 
I 4- 0(&). For better comparison the GUE curve is given in 
the re-scaled coordinates so that the asymptotic behavior at 
small and large k coincides with that of the GOE curve. 

More interesting is the asymptotic behavior of G(fc) at 
large |k| ^ 1. While Re[G'(K)] tends to a constant 



in Fig. |4] and Fig.jsjthe real and imaginary parts of G(k) 
together with the asymptotic behavior Eq. ([t]) at a small 
d2 = 0.01. The choice of only one constant, c = 1 + O(fe), 
ensures both the right leading logarithmic behavior of 
Im[G(K;)] and the correct limit of Re[G(K)] as S — oo. 
It is also seen that upon a proper re-scaling of x and y 
coordinates the shape of G'^°^(k) and G'^^^(k) is very 
similar. We conclude therefore that the shape of the G(k) 
dependence is dictated mostly by criticality encoded in 
the Chalker's conjecture and not by the symmetry of the 
Hamiltonian. 




FIG. 5. (color online) Imaginary part of G(fc) for GOE (blue) 
and GUE (red). The latter curve is given in the re-scaled 
coordinates as in Fig. [4] The black dashed line represents 
the asymptotic behavior Eq. ([7| for the proper choice of the 
coefficient c = 1 + 0.90 6 and = 0.01. 



Re[G(S)] = <( 2l{2' ^ 1 
/3 = 2 



ILL 

2 



|Im[G(K)]| increases logarithmically 



\A » 1), 



(44) 



VII. LEVY FLIGHTS FOR STRONG 
MULTIFRACTALITY: A POOR-MEN 
DERIVATION. 



Im[G(K)] 



/•It: 



ln(|K|) -I- const, /3 = 1 
-\ Ind^l) + const, /3 = 2 



(45) 



Such a behavior is fully consistent with Eq. ([t]) , in which 
a sub-leading term is added: 



/(x) = cx' 



(46) 



For very large x » exp(l/6), Eq. (461 reduces to Eq. (7) 
while in the intermediate region \ <^ x <^ exp(l/6), 
which is present only at small 6^1, one obtains: 

j{x) ~ cd2 log(a;) . 

Thus we confirm the Chalker's conjecture for the den- 
sity correlation function to the first order in the small 
parameter 6 <C 1. 

In order to illustrate the similarity of the density corre- 
lation functions conjectured in Eqs. ([5]), ([7| and (11 1 and 
obtained analytically for the CPLBRM Eq. (fT4|, we plot 



In this section, we exploit another expression for the den- 
sity correlation function (at a fixed energy of the wave 
packet £ which has been reinstated in this Section) in 
the time domain which is equivalent to Eq. ([T]) provided 
that w^i?orf;S>l/i?, where B is the energy band- 
width: 

V„,(t) = 2nY,{S{e - £„)e"(='"-^") (47) 

Note that Vr^r' (t = 0) = due to completeness of the set 
of normalized wave functions. Therefore, one can replace 
git(e„-e™) |-,y git(e„-em) _ ^ without changing the result. 

Let us consider very strong multifractality and very large 



distances R = r — 
approximated as 



Then the wave function can be 



(48) 
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where Hnm is an off-diagonal entry of the matrix Hamil- 
tonian and e„ = Hnn is the diagonal entry. In the case 
of very small b, when the wave functions normalization 
is almost all concentrated at one site n (center of local- 
ization), the eigenenergy of a state e„ is almost equal to 
the on-site energy Hnn- 



The main contribution to Eq. ( 47 1 is given by two terms 
(m = r, n = r + R) OT {n = m = r + R) . In both cases 
the combination of wave functions is equal to: 



\H, 



r,r+R\ 



Now, we average over disorder. For a Gaussian ensemble, 
averaging over £„ = Hnn and over the hopping matrix 
elements Hnm are independent from each other. There- 
fore, we replace |i?„mp by its average {\Hnm\^) given by 
Eq. ( |14[ ) and reduce averaging over e„ to the energy in- 
tegral with the help of the spectral correlation function 
R{e,e') = {p{e)p{e')): 



V. 



r,r-\-R 



it) 



de'Rie,e') . (49) 



Note that due to level repulsion R{e, e) = the integral 
in Eq. (49) is convergent at e' = e. However, at t being 
smaller than the Heisenberg time, t <C = Npo (or 
w ^ A ~ B/N), the region of level repulsion is very 
narrow so that one can replace R{e, e') « p{e)p{e'), where 
p{e) is the averaged DoS. Using also the symmetric form 
of the averaged DoS p{e) — p{—e) we obtain at e = 0: 

, , 47r/9o h'^ f^°° , , , ,x 1 - cos(te') , , 
-Dnr+Rit) « — 1^ de' p{e') (50) 



At t ^ the integral in Eq. (50) is dominated by 

e' <^ B and one can replace p(e') ~ po and obtain: 



cos{x) 27r2pQ &2 



i?2 



(51) 

This is exactly the result Eq. ( 26 ) obtained in the corre- 



sponding limit R^ ht from general formulae of the virial 
expansion. 



Note, however, that Eq. (51) does not describe the be- 



havior of the correlation function Eq. (47) at smallest 
t ^ Po ^ 1/B. In this case, one can expand the expo- 
nent in Eq. (50) and arrive at: 



i?2 



+ 00 



dxp{x) 



i?2 ■ 



(52) 



The two asymptotics Eq. (51) and Eq. ( |52| match each 

other dX t ^ pn. ^They both describe the Levy flights. 

However, Eq. (51) has the scaling form R~'^{t/R^'^) 
(at d — 1) while Eq. (52) is the fully perturbative re- 
sult which violates the critical scaling law 'Dr,r+R{t) — 
R^"^ S{R'^/t). In the macroscopic thermodynamic limit, 
when R — )■ oo, t/R"^ fixed but R/L — 0, the perturbative 
region plotted as a function t/R"^ shrinks to zero. 



VIII. DENSITY CORRELATION FUNCTION IN 
THE WIGNER-DYSON LIMIT b ^ oo 

One can consider another limit when the size of the sys- 
tem L = is fixed and the parameter b in Eq. (14) is 



increasing. The limit 6 — !• oo corresponds to the Wigner- 
Dyson RMT. In this limit, the system becomes effectively 
zero dimensional and is not critical any longer. Neverthe- 
less, we derive and briefly review its density correlation 
function for the sake of completeness. 



We start with the same Eq. ( 47 ) as in the previous section 



but employ the independent averaging over eigenvalues 
and eigenfunctions. The former is given by the famous 
Wigner-Dyson statistician -^^hile the latter is described 
by the Porter-Thomas statistics. The simplest case is 
/3 = 2 when the eigenfunction ipnif) = Unr is uniformly 
distributed over the unitary group U{N). The averages 
of the products of -^-functions are well known, e.g.: 

{UnriUn'r^^mr'J-Jm'r'-) = ^2 _ ^ ^^'^^ 
{^7vn' ^rir2^mm' ^r'-^r2 ^nm' ^rir^ ^mn' ^r2r2 ) 



Using Eq. (53), we obtain 



i27rpo)-^Vr,r+R{t) = S R,o + iK (t) - N) 
where K(t) is the spectral form-factor. 



Sr,oN ~ 1 
- 1 



, (54) 



Npo ^ 



At t = 0, summations over n and m are independent and 
one finds K[t = 0) = N while at large t oo only terms 
with n = m survive and K{t — > oo) = 1. However, the 
behavior of K{t) at small times is highly uneven: K{t) 
drops to almost zero for a very short time of the order 
of the inverse bandwidth po ^ 1/B oc 1/VN and then it 
recovers to unity; for example, in the unitary ensembli^ll 
K = t/tn for pq <^ t < ta and X = 1 for t > tn; we 
remind that the Heisenberg time is t-a = Npo. 
One can immediately see that at t = Eq. 
discrete (5-function Sr^q and that the sum rule Eq. 
fulfilled at any t. 
function remains 



( |54[ ) is the 
(g is 

Furthermore, the density correlation 
a combination of a (5-function and a 
flat background at any t. The weight of the ^-function 
determines the return probability: 



Pit) = 1 



N 



- 1 



(N-Kit)). 



(56) 



The non-monotonic behavior of K(t) results in the sim- 
ilar behavior of the return probability: for a short time 
t ^ Pq ^ 1/ VN, the wave packet leaves the origin almost 
completely. However, at later times it accumulates again 
in the origin reaching the value P{t) = at t > t^i- 
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Such an echo behavior is typical for a chaotic quantum 
system of finite size. 

The flat background behaves with time as 



27rpo Vr,r+R{t) 



N- K{t) 
iV2-f ' 



{R>0). (57) 



It rapidly grows at small times, then decreases a bit and 
reaches the constant limit 1/{N + 1) at t > t^. 
Despite all unevenncss of the behavior discussed above, 
there is one rough feature: at finite R the flat background 
is always small like N~^. Given that the critical Hamil- 



tonian Eq ( 12 ) approaches the Wigner-Dyson RMT at 
6 — TV, one expects the density correlation function to 
decrease at constant (large) system size N as the param- 
eter 6 > 1 increases. 



IX. NUMERICAL RESULTS 

In order to check the analytical results of the previous 
sections and, more importantly, in order to verify their 
robustness for long-range critical random matrix models 
with weak multifractality, we performed statistical anal- 
ysis of eigenvalues and eigenf unctions obtained by direct 
diagonalization of large matrices drawn from the Gaus- 
sian orthogonal ensemble Eq. (14) with (3 — 1. The 
results are summarized in Fig. [6| and Fig. [7] 




FIG. 7. (color online) Scaling function S — r £'r',r'+r(i) av- 
eraged over r' for different fixed r for the long-range critical 
random matrix ensemble Eq. (14 1 with /3 = 1, b = 5 and 



the matrix size iV = 10". Main plot: r = 200,500,1000. 



The region of Levy flights with scaling Eq. 1 13 1 shrinks com- 
pared to the case of small b = 0.1. The region of perturbative 



Levy flights Eq. (521 increases with decreasing r. Inset: the 
same scaling function for small r — 5, 10, 20 and larger values 
of t/r. The multifractal regime, Eq. ([8|, is well seen and it 



crosses over directly to the perturbative Levy flights, Eq. ( 52 I, 
at small t/r 



1/2 2 ^ 

b = 0.1, d=2 b + 0{b ); 


(t/r) = 


: N = ld' 






^T'^'t/r 








« r= 100 






" r = 200 






— r = 500 
" r= 1000 






,1 





t/r 



FIG. 6. (color online) Scaling function S = r I'r',r'+r(i) av- 
eraged over r' for different fixed r = 100, 200, 500, 1000 for 
the long-range critical random matrix ensemble Eq. ( |14[ ) with 
/3 = 1, b — 0.1 and the matrix size A*' — lO'*. The curves 
for the multifractal regime Eq. (|8| and for the Levy flights 
Eq. ( 13 1 collapse to the one single scaling curve. The pertur- 



bative Levy flights Eq. ( 52 1 at the smallest t do not obey the 
scaling. 

In Fig. [6| we demonstrate the scaling Eq. (|6| in the two 
main scaling regimes Eqs. ([8|, (13) for the case of strong 
multifractality b = 0.1. Both the multifractal regime of 
Chalker and the Levy flights oc t/r^ are well seen and 
the scaling Eq. ^ is confirmed by a collapse of curves 



for several values of the distances r to the one single scal- 
ing curve. The scaling Eq. (|6| is violated at very small 
t < pq where inhomogeneity of spectrum (finite size ef- 
fects in the energy space) should be taken into account 
and Eq. ( 52 ) stands for Eq. ( 13 1. Another possible source 



of violation of scaling is the finite-size effects in the coor- 
dinate space which exhibit themselves as the corrections 
in the parameter r'^/N ^ 1 and are important at not 
very small r'^/N. 

In Fig. [T] we show the critical density correlation func- 
tion for weak multifractality in the case 6 > 1. Besides 
the scaling and the Chalker's ansatz for the power-law 
behavior Eq. ([8| in the multifractal region, the focus of 
our study was the power-law tail at small t/r that de- 
scribes the Levy fiights. It is seen that the region of the 
"scaling" Levy flights Eq. ( 13 1 is much more narrow for 



6 > 1 compared to the case of small b and is shrinking as 
r decreases. However, the exponent p of the power-law in 
the scaling function S{x) cx x^^, (x 3> 1) stays constant 
p — 1 and independent of the strength of multifractal- 
ity, i.e. of the parameter b. This result is not obvious 
for large b where the arguments behind the derivation of 
Eq. (51 1 no longer apply. It was also verified analytically 
by considering the next order in the virial expansiorP^. 
In the next section, we will consider the classical analogue 
of such a behavior. 
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X. RANDOM WALKS OVER HIERARCHICAL 
MANIFOLD WITH LEVEL-DEPENDENT 
ASYMMETRIC RATE 

In Ref.|18j Wegner and Grossmann (WG) suggested a 
classical random walk process on an hierarchical struc- 
ture that is rich enough to mimic many relevant regimes 
in transport in disordered systems and in fluid turbu- 
lence. Its rigorous definition is given in the original work. 
Here we illustrate this process for the 2d Sierpinski gas- 
ket, see Fig. [8j As any regular fractal, the Sierpinski 
gasket is characterized by the hierarchy of self-similar 
clusters, each one containing z clusters of the next gen- 
eration (level). For the Sierpinski gasket of Fig. [sjz = 3. 
Another important parameter of the fractal geometry is 
the space scaling factor /i > 1 that is the ratio of linear 
sizes of the clusters of the fc-th and the fc + 1-th levels. 
For the Sierpinski gasket of Fig. |8]/^ = 2. As in the orig- 
inal work ^18) . we assume that there is a smallest level 
fc = with the size of triangle equal to 1. Then, the 
size of the triangle of the largest level k = I is L ~ 
and there are white triangles of the smallest level in it. 
The parameters z and /i determine the Hausdorff fractal 
dimension, i.e. the exponent dh that governs the scaling 
S cx L'''' of the total area of white triangles as the size of 
the largest triangle L increases. Given that S (x z\ one 
immediately finds that 



In z 
In n 



(58) 




FIG. 8. (color online) The example of a 2d hierarchical man- 
ifold: jumps occur only between the full (colored) triangles 
of the neighboring generations (levels) shown by different col- 
ors. There are neither jumps between triangles of the same 
level (size), nor jumps between triangles that are not directly 
touching each other. 

The WG random walk process can be viewed as jumps 



over the manifolds of centers of the full (colored) triangles 
each of them being determined by two numbers (fc, i), 
where k is the level of the triangle and i enumerates in a 
certain way (see Fig.jsJ all the z*^"' triangles of the given 
level k. An important rule is that a jump may occur only 
by one level up and by one level down. The up-jump may 
occur from each of z smaller full triangles of the k — 1 
level to the one single larger triangle of the A:-th level 
which smaller triangles touch directly. It happens with 
the rate: 



(59) 



The opposite process of down-jump is supposed to have 
a rate: 



W, 



(down) 



S 

z r 



(60) 



The presence of z in the denominator is natural as the 
down-jump happens to one of the z smaller triangles that 
touch the larger triangle. What is the crucial invention of 
Wegner and Grossmann is the rate asymmetry parameter 
r. It is this parameter that makes the up-jumps (down- 
jumps) a predominant process if r > 1(< 1). 
The rate scaling parameter s is natural on every hierar- 
chical structure of the type Fig.[8]in order to compensate 
for the increasing length of jumps as the level increases. 
The "normal" case is when s < 1. In this case the Levy 
flights cannot occur. On the other hand, if s > 1 and 
r — 1, all the moments {x"^) {m > 0) of random walk 
displacement are divergent at all times. This is a some- 
what pathological model. 

The presence of the rate asymmetry parameter allows 
one to reach new regime (regime "C" in Ref.jTnj) at ?■ < 1 
where s > 1 but (r s) < 1. It is exactly the regime where 
the density correlation function behaves similarly to that 
in the quantum case of random critical long-range Hamil- 
tonians of Eq. (12). Note that the case r < 1 and s = 1 



would result in trapping of a random walker on a frac- 
tal with no restriction to the level number from below. 
Indeed, in this case the bias towards down-jumps would 
drive the walker to lower and lower levels of hierarchy 
without a considerable displacement in the space. It is 
not the case in the WG model where the minimal level 
of hierarchy does exist, so after reaching it, the walker 
must go to the higher level with certainty and thereby 
may escape to arbitrary large distance. Thus, the coarse- 
grained WG model never leads to localization but only to 
sub-diffusion: 



(e<i). 



(61) 



In Eq. (61 1 we introduced the walk dimension d^ which 



is essentially the dynamical exponent that determines the 
relative scaling of space and time. In the WG model with 

(62) 



s > 1, 



sr < 1, 



it is given b 



In(sr) 
In/x 



(63) 
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and leads to a sub-diffusion < 1 if (sr) 



-1 >m'. 



where dh and d,„ are the Hausdorff and the walk dimen- 



The definition Eq(61) is similar to the celebrated critical 
subdiffusion (x)^ oc t^/"^ which is well-known from the 
scaling theory of Anderson localizatiorP^. This suggests 
that d-uj in (61 1 has the same meaning as the dimension- 
ality of space in the quantum critical models. Now we 
show that the WG random walk with the parameters s 



and r obeying Eq. ( 62 ) reproduces all the regimes shown 
in Fig.[l} 

First of all, we note that, according to the Table 3 of 
Ref. [in], the moments (|a:|'") are divergent if ^'"(rs) > 1 
which requires m > d^ at rs < 1 < s; here d^ is given 



by Eq. ( 63 1 and x is the x-coordinate of the displacement 
vector R = {x, xi, X2, ...Xd-i} in the d-dimensional 
space in which the fractal is embedded. This suggests 
that the probability to find a walker at in an interval 
(a;,x -|- dx) at a time t is proportional to dx. 
Switching to the radial-angle variables we find that the 
probability to find the walker at a distance interval 
{R, R + dR} within the solid angle interval dft is: 

P{R, t) cx R'^-^dRdn (64) 

(X R-"^ R-"^^- d'% R^oo. 

In order to find the coefficient of proportionality in 
Eq.([64| we use Eq.(4.30) of Ref. [121 . Then we obtain 

^^{rs)'t^ (65) 



R'^P{R,t) = S{R,t) 



1 



Here k = (1 — r)(l — rs) and I = 1{R} — In R/ In fi is the 
displacement in the ultra-metric space of levels expressed 
in terms of the radial displacement in the real space. At 
large enough R and sr < 1 the exponent can be expanded 
and we obtain: 

R-'^P{R,t) = SiK^^" It) = Kti?-'^?^^ = KtR-'^^. 

(66) 

We see that, in full agreement with the meaning of d^, as 
an exponent of dynamical scaling, R~'^"' enters the func- 
tion S{R,t) in the scale-invariant combination R'^'^' /t. 
Thus Eq. ( 66 1 , which describes the Levy flights for the 



classical WG random walk, is isomorphic to Eq.(13) for 



the density correlation function in the quantum prob- 
lem of random critical long-range Hamiltonians, provided 
that: 

d^ d. (67) 

Now consider small distances R'^"' ^ t. At such distances 
Eq.( |66| ) no longer applies. In order to find the density 
correlation function in this regime we apply the scaling 
S{R,t) — S{R'''^" /t) and the coarse-graining. The latter 
implies that the return probability P{t) is proportional 
to 5(1, i) = S{l/t). So, the time-dependence of return 
probability allows (via scaling) to find the entire function: 

SiR,t)^S{R'^^/t)(xP{t/R'^"'). (68) 

The return probability in the WG random walk process, 
regime Eq.(|62|), has been found in Ref. |20) : 

In z dh 
ln(rs)-i 



p{t) cx r", f 



(69) 



sions give n by E q. ( 58 1 and Eq. ( 63|) 
and Eqs.( [68l69l ), one finds for 
function at i?""" ^ t: 



Then using Eq.(64|, 



le density correlation 



R-'^S{R,t) cx R 



-d 



(70) 



This is a classical counterpart of Eq.Q. Quite naturally, 
the correlation fractal dimension d2 is replaced by the 
Hausdorff dimension of the classical fractal: 



dh -5> d2, 



(71) 



and again the correspondence Eq.(67) holds true. 
We conclude this section by saying that there is a com- 
plete quantum-to-classical analogy in the density correla- 
tion function of the WG random walks and the long-range 
critical random Hamiltonians. 



XI. CONCLUSION 

The main results of this paper are the following: 
(i) We have identified two regions with qualitatively dif- 
ferent behavior of the density correlation function of long- 
range critical Hamiltonians Eq. (12): the multifractal 
region where the power-law behavior Eq. ([s]) predicted 
by Chalkei^ is valid and the Levy flights region with 
the power-law behavior Eq. (13 1. Both types of behav- 



ior were studied analytically within the virial expansion 
method and by a direct numerical diagonalization of large 
matrices for the critical random- matrix model Eq. (14). 



(ii) It appears that for strong multifractality (6 <C 1 in 
Eq. ( 14 )) there is a complete analogy of the density corre- 



lation function in the quantum problem due to emergent 
fractality of random critical wave functions and in the 
classical random walks on hierarchical structures due to 
their geometrical fractality. In both cases, one can find 
two independent critical exponents in the density corre- 
lation function: Classical random walks on fractals can 
be described by the Hausdorff dimension dh and the walk 
dimension d^ (or the spectral dimension dg = 2dh/dw)- 
Two corresponding critical exponents of the quantum 
problem are the correlation multifractal dimension d2 and 
the dimensionality of space d. 

(iii) It is remarkable that only the latter determines the 
power-law large-distance tail (quantum Lev y fligh ts) in 
the scaling function S{t < r"*) cx t/r'^, Eqs.( |6|l3[ ). The 
exponent —d is independent of b (and hence on d2{b)) at 
all values of the parameter b in Eq.(14). By comparing 
Eq.([l3l) and Eq.(66), we find that the scaling variable 



X = r jt (quantum critical problem) or a; = r jt (clas- 
sical random walks) enters in the Levy flights power-law 
tail in both cases as x^^ with the universal exponent 
— 1 which is independent of the correlation dimension d2 
(quantum critical problem) or of the Hausdorff dimen- 
sion dh (classical random walks). We believe that this 
universality has a deep physical origin. 
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e-o L = 


10, 


r = 3"^ 


Q-^l L = 


20, 


r = 2x3'" 


0-^ L = 


30, 


r = 3x3'" 



t/r 



FIG. 9. (color online) Scaling function S — r^Vrfiit) vs. 
scaling variable t/r^ for the 3D Anderson localization model 
with the critical disorder W = 16.5 and the nearest-neighbor 
hopping integral V — 1. The plots for three different values 



of r = 3 



1/2 



r = 2 X 3 



1/2 



and 



= 3 X 3^''^, and the cor- 
responding three sample sizes L = 10, L = 20 and L = 30 
collapse (at t > tmax) in the one single scaling curve with the 
"multifractal" power-law behavior (r^/t)"^^/^, where d2 ~ 1-3 
is the correlation fractal dimension. The dotted-dashed line 
represents the corresponding power-law. At t < tmax there is 
a significant divergence of curves for different r. Given that 
the ratio r/L « 0.17 is constant for all the curves, this diver- 
gence can only be explained by a significant role played by a 
"microscopic length" (see Ref. [T5]) and the corresponding 
dimensionless combination {Lt/£o)^ = tEo ~ tV = t. Note 
that at maximum of the scaling curve t — tmax = 1.5, 12 and 
40 for the three values of r indicated above. For the same 
values of r, the smallest value of the scaling variable shown 
on the plot corresponds to t = 0.25, 2.0 and 6.8. The "macro- 
scopic limit" Lt/£o 3> 1 implies t 3> 1. Only in this limit one 
expects a collapse of the curves. This condition is obviously 
violated for the curve with the smallest r = 3^^^ at t < tmax- 



The application of the above results to the critical point 
of the Anderson localization transition is a subtle issue 
due to the absence of an exact solution to the Anderson 
localization problem and due to the limited sizes of 3D 
systems amendable to direct numerical diagonalization. 
Yet, it is rather well established that the multifractal 
behavior of the density correlation function, Eq. ([S]), is 
present at the Anderson transition point in 3D systems. 
In Fig. [9]we present the results of a direct numerical diag- 
onalization of the 3D Anderson model at critical disorder 
W = 16.5 and the nearest-neighbor hopping V — I. The 
density correlation function 2?r,o(i) is computed for three 
different values of r and the corresponding three different 
system sizes such that r/L « 0.17 — const. Thus geo- 
metrically the systems are "macroscopically similar" and 



the finite-size corrections of the type r/L, though appre- 
ciable, are the same for all the three cases. One can see 
that for t > tjnax ( tmax corresponds to the maximum of 
the scaling curve) all the three curves for r^ 'Drfi {t) col- 
lapse to the one single scaling curve S{r^ /t) ~ {r^ /t)'^^ 
with rather good accuracy. In this region they exhibit 
the power-law behavior which is represented by Eq. (|8|, 
albeit with the power slightly modified by the finite-size 
r/L effects. 

The situation with the tail at r^ ^ t is much worse. The 
curves for different r diverge significantly as the scaling 
variable t/r^ gets smaller. There is only one reason for 
such a behavior: it is the finite-size effects of the type 
io/Lt, where io ^ {poV)~^^^ is the microscopic length 
which has been introduced in Ref.|12|. and L^ = t/pQ. 
Here po ^ (a^ V W)~^ is the density of states, thus £o ^ 
a W^l'^ is of the order of the lattice constant a. The 
parameter that controls the smallness of such finite-size 
effects is {l^/Ltf ~ itV)-^ = l/t. Obviously, at small t 
this parameter is no longer small and one cannot expect a 
scaling behavior in the region t <^1. On the other hand, 
tmax ~ 0.3 r^. So, it is only at large r that one may 
observe a scaling behavior in a sufficiently wide interval 
1 ^ t < tmax- We conclude that in order to meet the 
requirement 1 <^t < tmax and simultaneously to exclude 
the finite-size r/L effects one should study really large 
3D systems with L ~ 100 which is numerically a very 
hard problem. 

In this situation we can only guess what is the form of 
the scaling function S{r^/t) at t < t„iax the Ander- 
son transition point of 3D systems. The most natural 
assumption is that it is exponentially small, and thus the 
Levy flights are absent, due to the short-range nature of 
the Anderson Hamiltonian. A counter-scenario in which 
the Levy flights may be present, is offerecF^ by the ob- 
servation of the " explosion of high-gradient operators" in 
the non-linear supersymmetric sigma-model in (2-f e) di- 
mensions. Unfortunately, this long-standing controversy 
still remains unresolved. 
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